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This work deals with thick braneworld models, in an environment where the Ricci scalar is changed 
to accommodate the addition of two extra terms, one depending on the Ricci scalar itself, and the 
other, which takes into account the trace of the energy-momentum tensor of the scalar held that 
sources the braneworld scenario. We suppose that the scalar held engenders standard kinematics, 
and we show explicitly that the gravity sector of this new braneworld scenario is linearly stable. We 
illustrate the general results investigating two distinct models, focusing on how the brane prohle is 
changed in the modihed theories. 
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I. INTRODUCTION 

This work deals with braneworld models in the pres¬ 
ence of scalar helds [in in an AdSs geometry with a 
single extra dimension of infinite extent. Differently from 
the original Randall and Sundrum (RS) proposal pre¬ 
sented in (H, which leads to the thin braneworld scenario, 
the coupling with scalar fields is also of interest and leads 
to another, thick braneworld scenario M- The pres¬ 
ence of a scalar held makes the warp function to behave 
smoothly, leading to a thick braneworld scenario. This 
possibility has opened a new area of study, and here we 
quote Ref. Q for some works on the subject. 

In the current work, we go beyond General Relativity. 
Before doing this, however, we should care about issues 
concerning the presence of degrees of freedom which may 
lead to instabilities in these theories, and here we rec¬ 
ommend the interesting reviews on the subject To 
go beyond General Relativity, we add two distinct con¬ 
tributions, one taking the metric tensor to change 
R = (the Ricci scalar, the trace of the Ricci ten¬ 

sor) to F{R), and the other, modifying Einstein’s equa¬ 
tion with the inclusion of auxiliary fields, as suggested re¬ 
cently in Ref. Q. To be more specific, we recall that the 
modification that changes R into F{R) in the Einstein- 
Hilbert action is older and was reviewed in Ref. Q. The 
modification introduced in the presence of auxiliary fields 
is more recent 0 and was studied in Refs. HO within 
the cosmological context, to see how the auxiliary fields 
may contribute to the cosmic evolution. Moreover, it has 
also been recently studied within the thick braneworld 
context, in five dimensions with a single extra dimension 
of infinite extent 0[i3- See also Refs. El El for other 
related studies on thick branes. 

In order to widen the interest in the work, here we 
add together both R and T = (the trace of the 

energy-momentum tensor) into F{R,T) to describe the 
gen eralized model, inspired by the recent investigations 
(la Ell. We then study how the two modifications act 
to change the standard braneworld scenario, in an AdS^ 
braneworld geometry with a single extra spatial dimen¬ 


sion of infinite extent. We implement the investigation 
in Sec. El where we discuss the general theory, including 
linear stability of the gravity sector. We then study two 
specific new models in Sec. uni working out the most 
important results of the models. As one knows, scalar 
fields play a fundamental role in cosmology as possible ex¬ 
planations for inflation, late time acceleration, and dark 
matter, among other issues. Also, scalar fields are fun¬ 
damental sources in the thick braneworld scenario with 
a single extra dimension of infinite extent. Thus, here 
we investigate scenarios described by a single real scalar 
field to source the thick brane configuration, leaving aside 
studies related to cosmology. We then end the work in 
Sec. lIVl where we include our comments and conclusions. 


II. THE PROBLEM 

We start investigating a model that describes general¬ 
ized gravity coupled to a scalar field in five-dimensional 
space-time, with a single extra dimension y of infinite 
extent. The model is controlled by an action of the form 

5 = I dy v^( - ^FiR, T)+Cs), (1) 

where R and T represents the Ricci scalar and the trace 
of the energy-momentum tensor, and Cs describes the 
source Lagrange density, which we suppose is of the form 

Cs = - U(<^). (2) 

It engenders standard kinematics, and we are using 
47rG^^^ = 1,5 = det{gab), a,b = 0,1,..., 4 and the signa¬ 
ture of the metric is (-1-). We focus attention on 

the prohle of the thick braneworld scenario, so the addi¬ 
tion of fermions and gauge helds will not be considered in 
this work. We take the spacetime coordinates and helds 
as dimensionless quantities, and we use Greek indices for 
the embedded (3-|-l)-dimensional space, y,!/ = 0, 1,2,3. 
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We can use the source Lagrange density © to write 
the energy-momentum tensor in the form 

"Tah = -^gabVc?5>V“(^-|-gab^ +(3) 

which contributes to give the trace 

r = -(3/2)V,</)V“</. + 5y. (4) 

Moreover, the equation of motion for the scalar field has 
the form 


A. Stability 


The investigation of linear stability of the braneworld 
model can be done assuming that the metric is perturbed 
in the form 


ds^ = 




h^j,v{y,x) dx^dx" - dy'^ 


( 11 ) 


Furthermore, the scalar field is written in the form 


(l) = (j){y)+^{y,x). (12) 


VaV> + + ^V^Ft + V^ = 0, (5) 


The first-order contributions in h and ^ that follow from 
the equations of motion lead us to 


where = dV/dcj), and Fr = dF/dR, Ft = dF/dT, 
etc. By varying the action of the gravitational field with 
respect to the metric tensor we obtain the modified Ein¬ 
stein equation 


FjlRab 2^9abF F ^^ab^ 


— VoVb iFr 


= 2Tab + -FtS/ a4‘'^b(j^ ■ 


( 6 ) 




+ ^e - df^d„h^ 

1 . 


Fr- 


-e ‘^^df,d^P{x, y) + ^h' F'^ = Q{x, y)y^u , (13) 


where h = h^'^. For simplicity, we have introduced the 
functions P{x,y) and Q(x,y), which are defined as 


We study the case of a flat brane, with the line element 

= e^^y^i.dxy'dx'^ — dy^ , (7) 

where A is the warp function, is the warp factor and 
rj^u is the 4-dimensional Minkowski metric, with signa¬ 
ture (-1-,—,—,—). We consider the case where both A 
and (f) are static and only depends on the extra dimen¬ 
sion, that is, A = A{y) and (j) = (j>{y). We then get 


( 



{4 + 3Ft)A' + ^F^]^' = (i + ^Ft)v^, 


( 8 ) 

where the prime denotes derivative with respect to the 
extra dimension. Also, the modified Einstein equation 
becomes 


P{x,y) = 


and 


h” + hA'h' + 




5E0C + 3(/)'^' + {(p'^ + 2V)h 


Frt , 


Frr 

(14) 


Q{x, y) = 3A'P'{x, y) + P"(x, y) - P{x, y) + 


12 


Fj' — 


^(A'h + + 

-Frt (4A'2+A")1 {(t>'^+2V) h 


( 4 A '2 + A")Frr - -Fr 


bA'h' 


(15) 


5 *“(i +jfr) = + 1.4'^ - jfS . (9a) 

As it is standard check, here we also note that the set of 
equations dH) leads us with the equation of motion (l8|). 

An important characteristic of the brane is its tension, 
which is given by 

r=Jdyp{y), (10) 


If we take F{R, T) —>■ R, in the case of General Relativity 
we get P{x, y) = 0 and 

Q{x,y) = 2V4 + 2^'e + lh" + ^A'h' 

+ (16) 

We can simplify the investigation of stability consider¬ 
ing the transverse traceless components of metric fluctu¬ 
ations, that is, we take 

= ^ = 0- (17) 


t) =2{A'\A")FR---2A'F'n.{9h) 


where p{y) = —e^^Cg is the energy density of the brane. 

The next important step is to check if the modification 
proposed above contribute to destabilize the geometric 
degrees of freedom of the braneworld model. We inves¬ 
tigate this issue studying linear stability of the gravity 
sector in the usual way. 


Thus, we can check that Eq. (IT^ reduces to the form 


(_ 52 _ 4 ^, 5 ^ ^ - ^dy) h,, = 


-2aFrT 


Fr L 






(5E^e+30r) ■ (18) 
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We see that the analysis of stability depends crucially on 
the function F{R,T), as expected. The problem simpli¬ 
fies if we assume that the function F{R,T) is separable 
in the form: F{R, T) = G{R) +F[{T). This choice allows 
to write the previous equation as 

(^-dl-4A'dy + e-^^a^^'>-^dy^h^, = 0. (19) 

We introduce the 2 -coordinate in order to make the 
metric conformally flat. We take dz = and we 

write 


z) = . ( 20 ) 


In this case, the 4-dimensional components of hy_i, obey 
the Klein-Gordon equation and the metric fluctuations of 
the brane solution lead to the Schroedinger-like equation 

{^■^ + U{z'^hy,i, =p^hy,v , ( 21 ) 


where 


U[z) = 


tA. 


3 djlnGji) 

2 ^ dz 


l/d{\nGR)\^ 1 d^Gt 


4 V 


dz ! 2 Gr dz'^ 
Fortunately, we can use this U{z) to write 


( 22 ) 


S' 1 'S' 




(23) 


where 


_d 3 1 d(lnGfl) 

dz 2 ^ 2 dz 


(24) 


The operator —d^ + U{z) is then non-negative, and so 
the gravity sector is linearly stable. This is an inter¬ 
esting result, which shows that despite the modifica¬ 
tion introduced with F{R,T), when we write F{R,T) = 
G{R) + F[ (T), the gravity sector of the braneworld is lin¬ 
early stable. We note that the stability behavior of the 
gravity sector only depends on the warp function, since 
R only depends on A{z) in this case. 


III. SPECIFIC MODELS 

Up to here, we used the generalized F{R,T) grav¬ 
ity coupled to a real scalar field to obtain the equa¬ 
tions of motion and show that the gravity sector of the 
thick braneworld scenario is linearly stable for F(i?, T) = 
G{R)+F[{T). The results motivate us to further explore 
specific models, investigating how modifications depend¬ 
ing on R and T change the standard scenario. We do 
this below, studying two distinct models. 


A. A simple model 

The first example we consider is described by the 
F{R,T) function 


F{R,T) = R-aT'^, (25) 

where a and n are real parameters. This is a simple 
choice, and now the two components of equation ([S]) be¬ 
come 


= -t ^ , (26a) 

A '2 = i 0'2 _ iy _ (26b) 

6 ^ 3 \ 12 ^ 

where T = + 5V. 

We concentrate on the simplest case, with n = 1. Here 
the two components of the Einstein equation becomes 




A f‘? d 3o( .n 4 5cr 

A'^ = - 6'^ - 

24 ^ 12 


V. 


(27a) 

(27b) 


To get to the first-order formalism, we introduce another 
function, W = W{(f)), which can be used to see the warp 
factor as a function of the scalar field. We do this writing 
the first-order equation 


A' = -1W, (28) 

where 7 is a real parameter to be determined. With this, 
we obtain that 




27 

4-3a 


W^. 


(29) 


Note that if 7 = 1 — 3a/4 the solution is identical to 
the standard case (a = 0). On the other hand, for a 
general 7 the solutions depend on the parameter a. Using 
equations ( 1 ^ and ( 1 ^ we can also obtain the potential 
as 


VU) = 

(4-5a)(4-3a) 3(4-5a) ' 


(30) 


Furthermore, the energy density can be written in the 
form 


p(y) = 


/lj2e2A(y) 4 _ 4 q , 




4—5a 

The energy of the thick brane becomes 


( 31 ) 


Fa 


ay^ 8 —9a 
(4-5a) (4-3a)2 


/ CO 

dy e^^(y^W^Wy )). (32) 

-00 


We note that the energy is controlled by a, and for 0 < 
a < 4/5 the energy of the brane is positive. 
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FIG. 1: The solution (I34II . depicted for a = b = 7 = 1 and 
oi — Q (solid line), a = 0.1,0.2, 0.3, 0.4 (dashed lines). 



FIG. 2: The warp factor depicted for a = b — ^ = 1 and 
a = 0 (solid line), and for a = 0.1, 0.2, 0.3, 0.4 (dashed lines). 


An interesting example which leads to analytical inves¬ 
tigation of its stability is given with the superpotential 
in the form 


W{(j)) = 2asin(b(/)), (33) 


where a and b are real parameters. With this choice, we 
obtain the solution ( 1 ^ as 


Hy) 


1 

— arcsin 
b 


tanh{By) , 


(34) 


where B = 4706^/(4 — 3Qf). In Fig. [T] we depict the 
behavior of the solution (l34|) . We note that the thickness 
of the solution is controlled essentially by a. Moreover, 
as expected, the solution goes to the asymptotic value 
(j) ^ (f> = ±7r/2b when y —>■ ± 00 . 

Using equation (E51) we can also get the warp function 
as 



-3;r/2 -n -njl 0 njl n 37 r /2 

FIG. 3: The scalar field potential, depicted for 0 = 6 = 7 = ! 
and a = 0 (solid line), q = 0.1, 0.2, 0.3, 0.4 (dashed lines). 


A{y) 


27 a 


In 


sech{By) 


(35) 


If a < 4/3 we can get the thin-brane limit when y —>■ ± 00 , 
in the form 


A{y —>■ ± 00 ) 



(36) 


Finally, the energy density (1^ can be written as 


p{y) = 


16 ( 70 ) 


2 r 


3(4-50)1 
sech(By) 


-l+l-f 


3B(l-a) 
70(4— 3 a) 


sech^(i3y) 


4’7'a 

■ 3 ^ 


(39) 


In Fig. [5] we depict the warp factor for some values 
of the parameters. We note that the warp factor narrows 
as a increases. 

Furthermore, the potential (l30l) can be written as 




(37) 


which is depicted in Fig. |31 Note that the potential has 
several global minima at tp = where 




16 ( 07)2 
3 (4 - 5a) ■ 


(38) 


Fig. |4] shows how the parameter a controls the energy 
density. The equation (l3^ can be integrated to give the 
energy of the brane in the form 


_ B /8-9aN F (s) 

“ “ “4^U-5a/ r(s + 1/2) 


(40) 


where s = 1 -|- 2'ya/3B. 

The behavior of the stability of the gravitational sector 
can be analyzed analytically for the particular case where 
B = 2'yal3. In this case we can make the change 

y = -^arcsinh(i32:). 

JD 


(41) 
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FIG. 4: The energy density, depicted for 0 = 6 = 7 = ! and 
a' = 0 (solid line), and for a = 0.1, 0.2, 0.3, 0.4 (dashed lines). 
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FIG. 5: The stability potential, depicted for B = 2/3 (solid 
line), and for 13 = 1,4/3,5/3 (dashed lines). 


With this we get the stability potential in the form 


U{z) = 


2-bB^z^ 
~ (l + B^z^f ■ 


(42) 


It is depicted in Fig. [5l for several values of the parame¬ 
ters. 

The above results show that the thick brane narrows as 
a increases, but the brane profile is qualitatively similar 
to the case of a standard brane. 


B. Another model 

Let us go on and investigate another model, this time 
with contribution from the Ricci scalar too. We assume 
that 


F{R,T) = R + PR^ -aT, (43) 

where a and /? are re al p arameters. The case with a = 0 
has been studied in [^, and now we add both the R^ 


and T contributions, to get analytical results. In Ref. 0 
one noted that the term R^ contributes to generate an 
interesting effect, the splitting of the brane, and here we 
want to see how it works. In this new model, the brane 
equations become 


4-3a 


0'2=a"+ ^ (2A""+ 16A"2+ 8kL4"'+ 5 A' 24 "), 


6 


(44a) 


4 3cr .9 d Scr ..n 
-V = A'^ + 


24 ^ 12 

+ M (^5A'4+ 3A'A"'+ 32A'^A"- 4A"2) . (44b) 


We note that if /3 = 0 we get to the equations ((27l) . Also, 
for /3 0 the equations involve the fourth derivative of 

the warp function, which makes it difficult to search for 
solutions. However, as suggested in if we choose the 
warp function A{y) as 


A{y) = In sech{ky) 


(45) 


where fc is a positive parameter, this allows us to write 
the relations (IT41) as 


V{y) 


4 — 3a 4 — 3a 
12fc2 ^ 15A:^ 

4—5a 4—5a 


(29 - 495'^) , (46a) 

1^(10-1455^-^1475^) , 

4-5aV /’ 


(46b) 


where 5 = sech{ky). We use equation (I46ap to see that 
j3 should be limited to the interval 


3 

232fc2 


</3< 


3 

160fc2 ■ 


(47) 


Moreover, we suppose that the parameter a is positive, 
smaller than 4/5. For this model, we can write the energy 
density as 


p{y) 


12 fc 2 12P(6-5a) 4 8/3fc^ ^ 

4—5a (4 —3a)(4—5a) 4—5a 

116(6-5a) ^2 I 98(8- 7 a) ^ 4 - 



(48) 


A closer examination of the energy density (H51) shows 
that it start to split when 


3(8- 7a) 

16(125- 116a)fc2 ■ 


(49) 


The splitting depends on a, so the modification intro¬ 
duced with the presence of T contributes to control this 
effect. In Fig. [B] we depict the energy density for A: = 1, 
a = 1 /2 and /? within the interval (j47|) . In this case the 
split starts approximately at /3 = 0.012. 

The solutions of equation (I46al) can be written as 


^(y) = ± 


6 


4-3a 


320^fc2 
4-3a 


EllipticE 


392^/c2 . 

160^/c2-3/ ’ 
(50) 
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FIG. 6: The energy density, depicted for k = 1, a = 1/2 and 
,9 = 0 (solid line), 9 = 0.002,0.006,0.010,0.014,0.018 (dashed 
lines). 
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FIG. 7: The stability potential, depicted for fc = 1 and j3 = Q 
(solid line), /3 = 0-002, 0.006, 0.010, 0.014, 0.018 (dashed lines). 


where (ps = arcsin 


tanh(fcy) 


and EllipticE is the elliptic 


integral of the second kind. For /3 = —3/232fc^ and /3 = 0 
we obtain the solution as, respectively. 


<^( 2 /) = 

/116-87a 


/ ® 

V 4 - 3a 

tanh(fcj/) 


(51a) 

(51b) 


We now study linear stability of the gravity sector. For 
this we make the following change of variables 

y = ^arcsinh(fc 2 ;). (52) 

k 

It allows to get the stability potential as 


where, f{z) = 1 + — 8/3fc^(2 — 5fc^z^). This potential 

is depicted in Fig. [71 The appearance of a lump in the 
well confirms the splitting of the brane as /3 increases. 
We note that the stability potential does not depend on 
a, although a contributes to split the brane, as shown 
above in (|4^ . for instance. 

The results show that the presence of controls the 
splitting of the brane, as also noted in Ref. ; here, 
however, the T term also contributes, adding quantitative 
modification to the braneworld profile. 

We can also consider adding a quadratic term in T 
with a quadratic term in the scalar curvature. In this 
case, we have 

F = R + - oT'^ . (54) 

The equations (jH) become 

-ir+arT = 

A" + ^p{2A""+lQA"^+SAA"'+bA'A"^ , (55a) 

V{P)-\p'^+^-^rT-jT^ = 

-3A'2 - 4,p{bA''^+ %A'A"'+ i2A'^A"- AA"^^ .(55b) 

To solve the model we assume that the warp function is 
yet given by gSl), and we take k = 1. This problem is 
much harder then the previous one, but we then rewrite 
the above equations in the simpler form 

= ^ 5-2 4f](29S^ - 49S''‘), (56a) 

-3 + 3S'2 - 4/3 (5 - 58S'2 + 495^). (56b) 

These equations have analytical solutions for appropri¬ 
ate choices of the parameters a and /3. We have examined 
some specific cases, and they are all similar. Thus, we 
take, for instance, a = 0.02304 and /3 = 0.00433; the 
solutions are 

piy) = ±parcsin[tanh(?/)], (57a) 

Vip)) = -2.31489 -b 3.00269 cos^ , (57b) 

where p = 1.19576. The energy density becomes 

p(y) =-2.31489 S'2-b 3.71762 5^ (58) 

It is depicted in Fig |8l showing standard profile. The 
absence of splitting is due to the fact that the analytical 
solution is obtained for a very small value of /3. 


U{z) 


3fc2 2-5k^z^ 56/3fc4(l - 3fc2z2) 

4 (1 -b kPz'^Y (1 -b kPz'^Y f{z) 

3136A:®z2y32 


IV. COMMENTS AND CONCLUSIONS 

In this work we studied some thick braneworld sce¬ 
narios in an AdS^ warped geometry with a single extra 


{i + k^z^Yf(zY ’ 
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FIG. 8: The energy density, depicted for the solution (1581) . 


dimension of infinite extent, with the gravity sector ex¬ 
tended to accommodate the F{R,T) modification, cou¬ 
pled to a single real scalar field with standard kinemat¬ 
ics. We used the equations of motion to show explicitly 
that the gravity sector is linearly stable when we consider 
F{R,T) = GiR) + H{T). 

We investigated two distinct models, the first one with 
F{R,T) = R — aT, and the other with F{R,T) = 
R -I- f3R‘^ — aT. In the first case, the contributions that 
appear from the extra term T induces quantitative mod¬ 
ifications in the thick brane profile, without changing its 
qualitative behavior. In the second case, the presence of 
R^ works to split the brane. This effect was identihed be¬ 


fore in [13 , and the modification introduced with a T also 
contributes to the splitting, but it adds no new effect in 
the braneworld profile. We then conclude that the mod¬ 
ification of the form R F{R, T) = G{R) + H(T) does 
not destabilize the gravity sector of the thick braneworld 
with an AdS^ warped geometry with a single extra di¬ 
mension of infinite extent. Moreover, the presence of the 
term H{T) = —aT contributes to change the braneworld 
profile quantitatively, although it adds no new qualitative 
effect to the brane. If we add the term H{T) = —aT^, 
it seems that the effect is similar, but here, however, one 
needs to explore a larger range of parameters. This task 
requires numerical investigation, an issue that is out of 
the scope of the work. 

The present study drives our attention to some issues 
of current interest, one concerning the investigation of 
how the quantitative modifications that we have identi¬ 
hed in this work act to localize fermion and gauge helds 
within the thick brane [l^, [l^. Another, very natural 
continuation is related to the hybrid brane recently pro¬ 
posed, in which one changes the potential of the source 
held to make the scalar held solution compact , to see 
how it adds in this new scenario. This last issue is now 
under consideration, and we hope to report on it in the 
near future. 
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